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Abstract 

In this work we study the propagations of normal frequency modes for quantum hydrodynamic 
(QHD) waves in the linear limit and introduce a new kind of instability in a double-degenerate 
plasma. Three different regimes, namely, low, intermediate and high magnetic field strengths are 
considered which span the applicability of the work to a wide variety of environments. Distinct 
behavior is observed for different regimes, for instance, in the laboratory-scale field regime no 
frequency-mode instability occurs unlike those of intermediate and high magnetic-field strength 
regimes. It is also found that the instability of this kind is due to the heavy- fermions which 
appear below a critical effective-mass parameter (/i cr = \/3) and that the responses of the two 
(lower and upper frequency) modes to fractional effective-mass change in different effective-mass 
parameter ranges (below and above the critical value) are quite opposite to each other. It is shown 
that, the heavy- fermion instability due to extremely high magnetic field such as that encountered 
for a neutron-star crust can lead to confinement of stable propagations in both lower and upper 
frequency modes to the magnetic poles. Current study can have important implications for linear 
wave dynamics in both laboratory and astrophysical environments possessing high magnetic fields. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. INTRODUCTION 



Since the early developments jl-3|, more investigations have appeared on the problems 
of wave dynamics in quantum plasmas j^-lO], a degenerate dense and cold ionized matter. 
A quantum plasma is characterized by the state of plasma in which the thermal energy 
of the constituent particles fall below a well-defined energy called the Fermi-energy and 
the particle waves overlap causing a new type collective phenomenon through the Pauli 
exclusion rule 11] . It has been shown that, in a quantum plasma the nonlinear response 



;o plasma excitations may differ fundamentally due to effects such as electron nonlocality 
12 1, degeneracy 13], and spin-orbit magnetization [ijj]. The nonlocality effects due to the 



Bo 



im- force has been shown to lead to some criticality in the quantum hydrodynamics (QHD) 
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16] . It has also been revealed that the magnetohydrodynamics MHD waves substantially 



altered when the normal degeneracy is changed to the Chandrasekhar relativistic degeneracy 
\v\ in a Thomas-Fermi model and different degeneracy-regimes in a quantum plasma can 
lead to different QHD-wave behavior when a stron g m agnetic field is involved [18]. On the 
other hand, the newly developed extended MHD 19M21] confirms that the plasma spin- 
polarization can lead to even more peculiar features in a quantum plasma. It has been 
shown that, the plasma magnetization which introduces a negative like pressure term can 
have significant effect on nonlinear density excitation in quantum plasmas 
the Jeans instability 23] • More recently a fully relativistic 
the spin effect has been proposed by F. A. Asenjo et al. [2 
Dirac theory. 

In the presence of strong magnetic field the quantum plasma can undergo a dou- 
ble quantum-pressure mechanisms {25], called the double-degeneracy analogous to double- 



22] and even alter 



rydrodynamic model including 
based on particle-antiparticle 



adiabatic 



26] process, due to the electron spin and orbital quantization. It has be shown 



that such stress-anisotropy can lead to a transverse plasma magnetic collapse in strongly 
magnetized plasmas j^] . Because of the transverse electron momentum quantization under 
the magnetic-field in addition to normal electron degeneracy new energy binding mechanisms 
may occur in a quantum plasma 28] . Two distinct magnetization mechanisms interplay deli- 
cately in a magnetization process, namely, the positive Pauli paramagnetism and the Landau 
negative diamagnetism [29!]. It has been shown that, the dominance of each of these mecha- 
nisms can lead to quite different nonlinear behavior in a quantum plasma {30 ]. According to 
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the value of the electron effective mass there are regimes in which the plasma becomes para- 
magnetic or diamagnetic 3jj. A vary large value of effective mass-ration (m e /m*) has been 
recently observed in Bismuth compounds [32]. On the other hand, for the superconductor 
leavy-fermion systems this value can be very small giving rise to the quantum criticality 



39[. In this investigation we show that such heavy-fermion quantum-plasma systems can 
develop a new kind of instability in MHD wave propagations. We consider different regimes 
of magnetic-field/density to show the relevance of our investigation to wide environments 
from laboratory to the astrophysical dense plasmas. The presentation of the article is as fol- 
lows. The QHD plasma model including the spin contribution is introduced in Sec. [TTJ The 
linear dispersion of spin-induced magnetoacoustic waves are derived in Sec. IIIII and the dis- 
persion and instabilities are discussed based on field/density regimes. Finally, a conclusions 
are drawn in Sec. IIVI 

II. DOUBLE-DEGENERATE QHD MODEL 

The full quantum magnetohydrodynamic fluid set of equations consists of 
continuity, momentum, Poisson's and Maxwell electromagnetic equations which 



has been derived in Ref. [33] for a general collisional quantum plasma. In the 
following we focus on the simplified model of uniform (in space and time) mag- 
netization which eliminates the role of electromagnetic interactions. Therefore, 
our simplified two-fluid quasineutral plasma model consists of contributions from 
the generalized pressure tensor, P, nonlocality (quantum diffraction) effect called 
the Bohma-force, Fq, and the Lorentz force, F L . We further assume a collision- 
less condition for the plasma due to the Pauli blocking mechanism which greatly 
reduces the fermion collision frequency in a quantum plasma compared to the 
ordinary ones. For a singly-ionized collisionless quantum plasma the quantum 
hydrodynamic fluid equation governing the electrostatic perturbations, reads as 

^f- + V • (n^etl^e) = 0, U(i >e ) = [u x (i if ,),U y (i iel ),U z (i je )\, 

m/^ = Te(E + u i)6 x B ) - ■ P i>e (^e) + ^V^g, (1) 
V • E = 47rp. 

Considering the collisionless quasineutral uniformly magnetized quantum plasma 
consisting of degenerate-electrons and classical nondegenerate-ions, we further 



simplify the model by neglecting the pressure due to ions compared to that of 
electron degeneracy pressure. The arbitrary strength uniform magnetic field, 
B , is assumed to be directed along the x axis. Therefore, the complete set 
of simplified QHD equations describing the dynamics of electrons and ions in 
magnetized degenerate plasma, neglecting the Bohm-force on classical ions (due 
to the much lower momentum of ions compared to that of electrons), can be 
written as 



% + V • (T^eUi,, 



0. 



5 = -^ + Wa ( Ui xx), 

£ + MUe X *)] = £ V0 

Acj) = 4ne(n e — rij) rj 0, 



-V • P e (n e ) + 



ti 2 



2m e rrii 



(2) 



where, u^, u c ( ifi ) = eB /m^c and P e (n e ) are the ion/electron velocity, mass, cy- 

clotron frequencies and the electron pressure tensor, respectively. In a collisional classical 
plasma when the characteristic plasma frequencies are lower than the collision 
frequencies the plasma remains isotropic in all directions. However, in the pres- 
ence of high magnetic field when the plasma species collisions are rare the plasma 
isotropy breaks down and the plasma may possess different characteristics par- 
allel and perpendicular to the preferred magnetic field direction. In such a case 
the pressure tensor can be written as 



P± + (P\\ - P±)b x b x (P|, - P ± )b x b y (P|, - P ± )b x b t 

(P|| - P ± )byb x P ± + (P|| - P ± )byby (P|| - P ± ) by b Z 

(P|| " P±)b Z b X (P|| - P ± )b Z by P ± + (P|, - P ± )b z b z 



(3) 



with b=B/ \B\ [26|]. In the presence of a uniform magnetic field the corresponding 
tensor can be expressed as diagonal with element only parallel and perpendicular 
to the ambient field as 



P|| 
Pi 
Pi 

which introduces a local anisotropy leading 
Chew-Goldberger-Low (CGL) plasma model 



(4) 

;o the so-called double-adiabatic 
34} ]. Some authors have recently 



investigated CGL anisotropy effect on the wave dynamics in the classical plas- 
mas [35|, |36j. On the other hand, when the when the plasma species become very 
close to each other and the particle-waves overlap new collective features emerge 
due to the spin interactions and a classical plasma may develops quantum hy- 
drodynamic feature {s 1 ^ . It has been shown that in the direction transverse to 
the magnetic field quantum plasma experiences additional negative-like (com- 
pared to electron degeneracy force) magnetization force, F m = B Vr, where, T 
is the quantum plasma magnetization vector [38] which can be considered the 
total transverse pressure P e ± as P e ± = P e u — TB [281 ] (e.g. see the footnote in the 
article). In the reduced quantum plasma model where the motion of electrons 
are confined to x-y plane (V = [dx,dy,0]), the plasma pressure can be anisotropic 
perpendicular and parallel to the existing uniform field analogous to the CGL 
model. However, in this case the anisotropy is not due to real pressure but is 
due to the plasma magnetization having a purely quantum mechanical origin. 

Sondheimer and Wilson [401 ] have investigated the diamagnetism due to the 
electron orbit and paramagnetism due to the electron spin in a zero-temperature 
ideal Fermi-Dirac gas for an arbitrary magnetic field strength along with oscilla- 
tory terms leading to the known de Haas- van Alphen effect {29]. It is well-known 
that the plasma magnetization can be due to both electron spin and orbital 
quantization. 

Apart from the oscillatory term the Pauli spin-magnetization per unit volume 
is given as Tp = (3fj IB n e /2kBTp e ) B (with Tp e being the electron Fermi-temperature 
related to the electron number density in quantum plasma being much larger 
than the electron kinetic temperature). On the other hand, a diamagnetic con- 
tribution to fully degenerate plasma is caused by the Landau orbital quantization 
which introduces a new term to plasma magnetization Tp = — l/3Tp exists (this 
expression does not include the electron effective-mass). This leads to the total 
magnetization of T = T P + T L = n e /j, 2 B B /Ep e (E Fe = k B T Fe = (h 2 /2m e )(3ir 2 n ) 2 ^ is 
the Fermi-energy) for a ordinary quantum plasma. For many compounds the 
effective electron mass may differ substantially from that of noninteracting elec- 
tron mass, leading to a generalized expression for magnetization per unit volume 
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of the ideal zero-temperature quantum plasma, as 



r 

v 



kBT Fe 
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T, 



h 2 



Fe 



2kBm\ 



■(37r 2 no 



,2/3 



(5) 



with the spin-orbit magnetic susceptibility defined as 



Xs 



1 dT 

vim 



H%D(E Fe 
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(6) 



where, \l = m e /m* e is the fractional effective electron- mass parameter with m* 
being the effective electron mass and D(E Fe ) = 3n e /2E Fe the density of electronic 
states (DoS) at the Fermi level. It is noticed that the susceptibility may be 
positive or negative depending on the electron effective mass ratio in Eq. (TBT) r 
hence, the plasma may behave paramagnetic or diamagnetic. It has been shown 
that, in some Bismuth compounds the fractional effective electron mass can be 
as high as \i ~ 1000 |32| or some semiconductors. On the other hand, some semi- 
conductors or superconductor heavy-fermion systems may possess the effective- 
mass ratio values as low as \i ~ 10 3 [39]. In many semiconductors the value of 



the parameter, contributes to a measurable Landau diamagnetic effects 11]. 
Therefore, defining the Zeeman parameter as, e = fisBo/E Fe , the magnetization 
pressure contributes to the below negative transverse force (in y-direction) in 
the normalized momentum equation. 



Bn 
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In n. 



(7) 



Now considering the parallel component of the nonrelativistic electron degen- 

r In 

eracy pressure as P e \\ = (2/5)E Fe non e (with no being the equilibrium number- 
density of the electrons), we are led to the normalized and reduced set of QHD 
equations as below 



d t n + d x (nu x ) + dJnu 



y 

d t u x + (u x d x + u y dy) u x + 8 x n^ - H*£ \± (*g + 



0. 



Inn 

ay 



dx ^/n \ dx 2 dy 2 

d t Uy + (U X d X + Uydy) Uy + OyU^ ~ ^ ^ 1 ~ "" 

^ dy \Jn ^ dx 2 dy 2 J Qu z 0, 

8 t U Z + [U X 8 X + Uydy) U z + QUy = 0. 



(8) 
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Note that we have dropped the (i,e) indices for simplicity and have used the following 
scalings for the normalization 

It E 

V -»■ — V, t -» — , n( i)6 ) n n(i ie ), u (i) e) -+ c s u (i) e), (f>-> — -(f), (9) 
Af ujpi e 

where, w P j = ^Ane 2 no/mi is the plasma ion frequency, Aj = c s /uj p i is the ion gyroradius, and 
c s = ■>/ Epe/rrii is the quantum ion sound-speed of plasma. Furthermore, we have neglected 
the term proportional to small fraction m e /mj and introduced the electron nonlocality H = 
(Hujpi/ Ep e ) ^Jvriil (2m*) , the normalized Zeeman e = [ibBq/ Ep e , and the normalized ion- 
cyclotron frequency Q = eB / (crriiOJpi) parameters. 

The reduced equation set may be linearized using a Fourier-transform using a plane wave 
G = Go + GiExp[i(k ■ r — cut)) with G = [n,u x ,u y ,u z ], Go = [1,0,0,0], and k = [fey, 
(k 2 + kj_ — 1) to give the first-order perturbation relations, as 

— urii + h\u x i + k^Uyi = 0, 

-uu xl = -^fcp! - l k \\ n ii ^ 
-uuyi = -^-k^rii - \kx_rix + ^(1 - \)k\nx + ju z i, 
-uu z i = -ju yl . 

Consequently, the dispersion relation, obtained from the above linear equation-set can be 
calculated as 

fcj(4 + 3H 2 )(u 2 - VI 2 ) + oo 2 [k\{3H 2 + 3e 2 (/i 2 - 3) + 4) - Q(oo 2 - ft 2 )] = 0, (11) 

with the two distinct upper (fast) and lower (slow) modes, u± respectively, given as 

A = (4 + 3# 2 )-3A; 2 e 2 (3- y u 2 ) + 6ft 2 . (12) 

It is evident that for A 2 < 2Ak 2 (A + 3H 2 )tt 2 the two modes become unstable in Eq. ffT2l . 
This condition sets a criteria for the stability of plasma waves based on four independent 
parameters, namely, the equilibrium number-density, no, the magnetic-field strength, Bq, 
the fractional effective-mass parameter, /x, and the angle of wave propagation with respect 
to the magnetic- field, 9 = arccos/qi = arcsinA;_i_. The plasma fractional parameters may 
be related to these parameters as; e ~ 1.587 x 10 6 i?o^o A* -1 ) H ~ 7200. 5r?,Q 1 ^ 6 /i -1 / 2 , 

— 1/2 

and Vl ~ 7. 277 B n . In the proceeding section we will evaluate the stability of waves 
based on these independent parameters in different quantum plasma circumstances, such as 
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± 



12 



A ± . A 2 - 24k 2 (4 + 3H 2 )Q 2 



magnetized laboratory and astrophysical degenerate plasmas. It is tempting to find the 
classical dispersion relation from Eq. 1111 by letting H — y (or equivalently by 
letting H — > 0). However, as has been pointed out by Haas et al. [ 1 121 ] one can 
not obtain the classical limit of the dispersion relation in this way, since in the 
classical limit, the equation of state of the quantum plasma which has been used 
in this model also vanishes through the Fermi-energy which is proportional to 
the plank constant. 



III. POLAR MODE-DISPERSION AND HEAVY-FERMION INSTABILITY 



Figure 1 depicts the low-field polar dispersion curves for laboratory degenerate plasmas, 
i.e., (uq ~ 10 21 /cm 3 and -Bo — 30T). Figures 1(a) and 1(b) compares the variation of dis- 
persion for lower and upper modes of QHD waves due to change in the strength of ambient 
magnetic field. It is observed that, the lower mode frequency increases due to increase in 
the magnetic field strength, particularly, in direction perpendicular to the magnetic field 
direction, while, the upper mode frequency it stays unchanged. Figures 1(c) and 1(d) re- 
marked that, as the fractional effective- mass parameter increases the lower mode frequency 
does not change for /i < /i cr = a/3, while, it decreases for large increase in the value of this 
parameter. It is noted that, the fields of order ~ 40T may be easily accessible in high- field 
laboratories. It is further revealed that both lower and upper mode frequencies are always 
stable in low-field degenerate environments. 

Figures 2 and 3 shows environment of astrophysical degenerate (no ~ 10 22 /cm 3 ) under 
medium-field strength (Bq ~ 10 4 ) which may be present in magnetosphere of a highly 



magnetized star such as pulsar or magnetar 



411 ] . It is observed from Figs. 2 and 3 that in 



medium-field limit the situation is quite different from that in low-field limit shown in Fig. 
1. The increase in the magnetic field strength in this case affects both lower (Fig. 2(a)) 
and upper (Fig. 3(a)) frequency modes. In fact the increase of the magnetic field strength 
increases the frequency for lower mode, while, it decreases the frequency for the upper one, 
leading to instability in direction perpendicular to the field for some higher field strength in 
both frequency modes. From Figs. 2(b) and 3(b) it is further remarked that the increase 
in the fermion density has quite opposite effects on lower/upper modes compared to the 
variation in magnetic field strength shown in Figs. 2(a) and 3(a). 



S 



On the other hand, Figs. 2(c) and 2(d) remark distinctive behavior of polar dispersion 
of lower mode in intermediate-field regime. It is observed that, for /x > /i cr (Fig. 2(c)) the 
mode instability unlike for /i < /i cr does not occur. The instability of this kind occurs for 
very small values of the effective-mass parameter and is due to presence of heavy-fermions in 
the degenerate magnetized plasma. However, both effective-mass parameter region confirm 
that, the frequency of this mode decreases as the value of /i departs from the critical value, 
fj, cr = y/3. Moreover, Figs. 3(c) and 3(d) indicate that the upper mode follows similar 
feature in opposite direction, i.e., distinct differences in effective- mass regimes follow but 
the frequency in this mode increases as the value of \i departs from that of the critical. 
However, the instability is about to set-in for heavy-fermion /j-values. 

Figure 4 depicts the instability growth-rate in polar view for the medium-field regime for 
lower and upper frequency modes corresponding to the heavy-fermion limit. For the lower 
mode (Fig. 4(a)) it is remarked that, the increase in the value of fractional effective-mass 
parameter decreases the instability growth-rate in direction perpendicular to the external 
magnetic field. However, for the case of the upper mode the instability growth is inclined 
towards the perpendicular direction to the field decreasing slightly in magnitude as the 
effective mass parameter value increases in the heavy-fermion regime. 

In fig. 5 we have evaluated the extreme regime for high electron number-density (n — 
10 28 /cm 3 ) under high magnetic-field (B ~ 10 8 T). This situation may be encountered for 
a strongly magnetic neutron-star (NS) crust. Figures 5(a) and 5(b) indicate the similar 
features as Figs. 2(a) and 3(a) for the intermediate field strength but highly unstable lower 
and upper modes in the high field limit. The frequency change in this regime, however, is 
observed to be mainly confined to the magnetic poles. The polar heavy-fermion instability 
growth-rates for the two mode in this case are also resemble that of intermediate-field shown 
in Figs. 4 but somehow pronounced. It is concluded, generally, that critically high-fields 
present in dense NS crusts cause the magnetosonic waves to be confined in the direction of 
the ambient strong field, namely, the magnetic poles. 

It is evident that, the instability presented in this work is due to the interplay 
between electron spin and orbital magnetization. In other words, the Pauli spin 
magnetization tends to destabilize the plasma, whereas, the Landau orbital mag- 
netization has the stabilization effect. Therefore, the instability arises when the 
Landau orbital magnetization or the plasma diamagnetism becomes vanishingly 
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small compared to the Pauli paramagnetism. 

IV. SUMMARY AND CONCLUSIONS 

We studied the propagation the linear propagation of normal modes of QHD waves in 
a double-degenerate plasma using standard linearizing method. Three different regimes, 
namely, low intermediate and magnetic field strengths was considered and found that in 
the laboratory-scale field regime no instability occurs unlike those of intermediate and high 
fields. We have found a new instability affecting both lower and upper mode frequencies 
appearing in the heavy-fermion (// < y/3) limit so that the responses of the two modes to 
fractional effect-mass change in this regime is quite opposite to each other. Furthermore, it 
was shown that, the heavy-fermion instability due to extremely high magnetic field such as 
that present in a neutron-star crust can lead to propagations confines to the magnetic poles. 
This study can be appropriate for investigating wave dynamics in the linear limit for both 
laboratory and astrophysical environments possessing high magnetic fields. 
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FIGURE CAPTIONS 



Figure- 1 

(Color online) Variation in polar dispersions for lower- and upper frequency modes for the 
laboratory scale low density and low magnetic field-strength regime with respect to change 
in various plasma parameters, namely, the magnetic field strength,!?, and the fractional 
effective-mass parameter, /i, with variation in one plasma parameter at each plot. The 
thicknesses of dispersion-curves vary according to change in the varied plasma parameter. 

Figure-2 

(Color online) Variation in polar dispersions for lower-frequency mode for the intermedi- 
ate scale magnetic field-strength and electron number-density regime with respect to change 
in various plasma parameters, namely, the electron number- density, n , fractional effective- 
mass parameter, //, and the magnetic field strength,!?, with variation in one plasma parame- 
ter at each plot. The thicknesses of dispersion-curves vary according to change in the varied 
plasma parameter. 

Figure-3 

(Color online) Variation in polar dispersions for upper-frequency mode for the inter- 
mediate scale magnetic field-strength and electron number- density regime with respect to 
change in various plasma parameters, namely, the electron number-density, no, fractional 
effective-mass parameter, /i, and the magnetic field strength,!?, with variation in one plasma 
parameter at each plot. The thicknesses of dispersion-curves vary according to change in 
the varied plasma parameter. 

Figure-4 

(Color online) Variation in polar instability growth-rate of lower- and upper-frequency 
mode for a intermediate scale magnetic field-strength and electron number-density regime 
with respect to change in the fractional effective-mass parameter, u. The thicknesses of 
dispersion-curves vary according to change in the varied plasma parameter. 

Figure-5 
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(Color online) Variation in polar dispersions and polar instability growth-rate for lower- 
and higher frequency modes for neutron-star crust scale magnetic field-strength and elec- 
tron number-density regime with respect to change in various plasma parameters, namely, 
the electron number-density, no, fractional effective-mass parameter, u, and the magnetic 
field strength,!?, with variation in one plasma parameter at each plot. The thicknesses of 
dispersion-curves vary according to change in the varied plasma parameter. 
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